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We find that, at order a s , the partial width of Z° to heavy flavors receives a power correction 
from a novel QCD mechanism, which is not suppressed by inverse powers of Mz, but only by two 
unknown 0(Aqcd/w) constants. The hadronic W width also receives a similar correction. These 
parameters may be fitted from the global electroweak analysis, and consequently the Standard Model 
predictions of various electroweak observables will be updated. This new mechanism is of no help 
to reconcile the discrepancy in b forward-backward asymmetry. We also point out the implication 
of this mechanism to heavy flavor production in other collider experiments. 



The precision electroweak measurements at LEP and 
SLC provide the most accurate knowledge of the fun- 
damental parameters in the Standard Model (SM), and 
meanwhile place stringent constraints on the new physics 
beyond SM [l|. Though the overall agreement between 
the measurements and the SM fits is acceptable, few 
acute discrepancies still persist 2]. For the Z-pole ob- 
servables, there are currently a 2.4 a discrepancy in 
the forward-backward asymmetry of b quark (A b FB ), 
and a 1.9 a deviation in the peak hadronic cross sec- 
tion (crjj ad ). While there was a significant deviation in 
Rb = r[6fe]/Fhad in past years, the latest measurement is 
consistent with the SM prediction at 1 a level. 

In order to search the possible hint of new physics, 
one needs first scrutinize all the uncertainties within SM. 
Ubiquitous (nonpcrturbative) QCD effects constitute a 
particularly important and difficult source of the theo- 
retical uncertainties. This is best illustrated by muon 
g — 2, where the current theory error is dominated by the 
low energy hadronic contribution ,3]. 

In this Letter, we will report a novel QCD mechanism 
which renders a new correction to the partial width of 
Z° to heavy flavors at order a s . It can be calculated 
in perturbative QCD, up to two unknown constants of 
0(Aqcd / m )- Incorporating this new effect will generally 
influence the current SM predictions of many electroweak 
observables, especially Thad, Oh ad , as well as Rb, R c - It 
will also affect the hadronic W width in a similar way. 
We also comment on the implication of this mechanism to 
the heavy flavor production in other collider experiments. 

The hadronic Z° width has been measured to per mille 
accuracy. To match such a precision, radiative correc- 
tions from both QCD and electroweak sectors must be 
computed in comparable orders. For instance, perturba- 
tive QCD corrections have been included to 3-loop order 
for massless quark, and b quark mass effects have been 
included up to order of mf/M^ |4|. 

Aside from the leading power contribution, the 
hadronic width of Z° are also affected by power cor- 
rections. For Z° decaying into light hadrons, the well- 
known power corrections are characterized by the quark 
and gluon condensates, whose effects are suppressed by 
a factor of 1/Mj. Therefore they can be neglected in 



accordance with present experimental precision. 

The situation for heavy flavors is dramatically differ- 
ent. Because the heavy quark mass sets a new scale, 
the most significant power correction may start at order 
Aqcd/tti, thus more important than those suppressed by 
powers of Aqcd /Mz and m/Mz- Surprisingly, this pos- 
sibility has been largely overlooked, perhaps due to the 
difficulty for the standard Operator Product Expansion 
(OPE) to tackle such a multi-scale problem. 

However, as we will see, the recently-developed heavy- 
quark recombination mechanism (HQR) [M H can 
realize such an 0(Aqcd/"^) correction to Z° hadronic 
width. In fact, when HQR was introduced and first 
applied to the B hadroproduction, it was shown that 
HQR generates an 0(KQc^/mb) correction to the total 
B cross section With a more concrete hadroniza- 
tion picture, HQR was originally motivated as a "higher 
twist" mechanism, to supplement the usual heavy quark 
fragmentation. However, by studying Z a decay into 
B at 0(a%), one recently realizes that under some cir- 
cumstances, HQR can overlap with the fragmentation 
mechanism ||. In that work, the 0(Aqcd/"^6) con- 
tribution from HQR must be identified with the con- 
tribution to the fragmentation function. Had we been 
able to extract the finite "higher twist" contribution by 
removing the "leading twist" term, it would represent 
an 0(o^AQCDWb/Mf ) ~ 10~ 6 correction to the partial 
width of Z° to bb, thus negligible with current exper- 
imental precision. In fact, the goal of this work is to 
show there is a new HQR process occurring at order a s 
only, with a genuine "higher twist" contribution of or- 
der Aqcd /rrib- Therefore, it is mandatory to consider its 
impact on the heavy flavor observables. 

The central picture of HQR is that, a heavy quark can 
bind with a light parton which emerges from the hard- 
scattering and carries an O(Aqcd) momentum in the rest 
frame of heavy quark. It is somewhat analogous to that 
a proton can capture an incident low energy electron to 
form a hydrogen atom. The word parton deserves some 
elaboration. Intuitively, one expects that b and q are 
more inclined to bind into a B meson, and bq diquark 
tends to evolve into a b baryon. Indeed, the bq y| and 
cq 0] recombination have been developed and applied 
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to a variety of fixed-target experiments to explain the 
observed charm meson and baryon asymmetries 0, Q ■ 

The last recombination mechanism awaiting explo- 
ration then is the bg recombination, when the parton is 
a gluon. At first sight, it seems in contradiction with the 
picture of constituent quark model. Indeed, this objec- 
tion is justified for the heavy quarkonium case, e.g., it is 
unlikely for bg to evolve into B c . However, the dynam- 
ics of heavy-light hadron is rather different from that of 
heavy quarkonium. Since the soft gluon in bg can eas- 
ily split into qq pair, and bg can easily pick up a light 
antiquark from the vacuum to hadronize, there is no dy- 
namical reason for bg recombination to be suppressed 
relative to bq and bq recombination. On the other hand, 
bg recombination does receive one parametric suppres- 
sion from the large N c consideration. There is a 1/N C 
suppression for bg to evolve into a B meson relative to 
bq, since it requires an extra qq pair. Similarly, it is even 
less probable for bg to evolve into a b baryon than B, due 
to a further penalty of 1/N C . 

At order a s , bottoms are produced through Z° — > bbg. 
This represents an ordinary 3-jet event if each parton in- 
dependently fragments. Nevertheless, in a small corner 
of phase space where g is soft in the rest frame of b, they 
can form an intermediate bg state with definite color and 
angular momentum. This state then hadronizes into a b 
hadron plus additional soft hadrons. Therefore, we end 
up with a jet containing a b hadron from the recombi- 
nation and a recoiling b. The corresponding Feynman 
diagrams are depicted in Fig. ^ We label the momenta 
of bg and b by p and p. 

The bq recombination respects a simple multiplicative 
factorization It is suggestive that the inclusive B 
production from bg recombination may be also written in 
a factorized form (The symbol B schematically represents 
any ground state hadron containing a b quark): 

T[B] = ]T t[Z° - bg(n) + b] i[bg{n) -> B] , (1) 

n 

where T n are the perturbatively calculable partonic cross 
sections, and are so-called recombination factors, 
which roughly amount to the probability for b and g to 
evolve into a state including B. The color and angu- 
lar momentum quantum numbers of bg are collectively 
labeled by n. The inclusive B production from bg recom- 
bination is identical to Eq. (|TJ, because of CP invariance. 

These £„ factors are analogous to those p n , rj n asso- 
ciated with bq and bq recombination. Thus far, a rig- 
orous definition in terms of nonperturbative matrix ele- 
ments is only available for p n yj, but generalization to 
rj n and should be straightforward. We also assume 
f£ ~ Aq C d/™&, the same as p n 

The bg that emerges from the hard-scattering can be 
in either of three irreducible color representations: 3, 6 or 
15. Fortunately, due to very simple color structure, only 




FIG. 1: Feynman diagrams for the bg recombination process 
Z — > bg(n) + b. Wavely, solid and curly lines stand for Z° 
boson, b and gluon, respectively. The shaded blob represents 
the hadronization of bg(n) into a B hadron plus anything. 
Note only diagram b) has a nonzero contribution. 

color-triplet channel survives in this process. The color- 
triplet state bg can be labeled by a fundamental SU(3) 
index i and has the normalized color wave function 

\bg(3)i) = ^TjlbJte*), (2) 

where SU(3) generators T a (a = 1 • • • 8) satisfy the nor- 
malization tr(T a Tfc) = S a b/2. In more complicated bg re- 
combination processes, the other two color channels will 
also contribute. 

We also need project the bg onto the states of def- 
inite angular momentum. We will consider only the 
S-wave states, because the higher orbital angular mo- 
mentum states will result in towers of Aqcd/w{, sup- 
pression. Therefore, the bg state can be either 2 Si/2 
or 4 S f 3/2, corresponding to a spin-^ or spin-| (Rarita- 
Schwinger) fermion. The momentum eigenstates of these 
fermions are conveniently represented by a vector-spinor 
u j(p) E3i which carries an extra Lorentz index relative 
to the Dirac spinor. Both of Uy 2 (p) and Uy 2 (p) satisfy 
the on-shell condition ftuj(p) = and transversality 
constraint p ■ uj(p) — 0. The spin-| vector-spinor is sub- 
ject to a further constraint 7^3/2 = ^" 

In addition to projecting the amplitude onto states of 
definite color and spin, we need isolate those most sin- 
gular terms when the gluon has small momentum in the 
rest frame of the b quark js|. This can be accomplished 
by setting p g —> x g p, where Xg IS cL momentum fraction 
that is presumably of order Aqcd/^&j and then taking 
the limit x g — > 0. The factor \/x g appearing in the am- 
plitude must be absorbed into the nonperturbative factor 
£„. Pragmatically, this can be described by a simple pre- 
scription. First we make the following substitution in the 
bbg amplitude: 

u 3 {p b ) e*^{p g ) - x g e[6 5 (3, J)]i ^ T a ) m b u»(p) , (3) 

then set p g = x g p in the rest of the amplitude and take 
the limit x g — > 0. 

To justify a perturbative treatment, we need show the 
parton process is really governed by the short distance 
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scale. The characteristic momentum carried by the re- 
combining gluon is p g ~ Aqcd^, where v = p/rrib is the 
4- velocity of bg. Therefore, the virtual b quark in Fig. Ok) 

travels a distance of order J A ^ P / Aqcd before turning 
on-shell, while the corresponding distance in Fig. ^d) is 
of order ^ "^f^ /AQCD, both are much smaller than 

the confinement length 1/Aqcd- 

Nevertheless, Fig.^i) is potentially worrisome, since it 
is reminiscent of a 5(1 — z) like fragmentation contribu- 
tion, of which we need get rid, to avoid double counting. 
The subamplitude of Fig. QJ,) involves a term 

Ujiph^ + m) = 2uj(p)-p-u^(p)(p-m)j tl . (4) 

Pleasingly, by transversality and the equation of motion, 
one immediately sees that it vanishes for both 2 S\j 2 and 
4 5 , 3 / 2 states. So only Fig.QJ)) needs to be included, which 
has nothing to do with the b fragmentation, since the re- 
combining gluon is emitted from b. Therefore it repre- 
sents a genuine "higher twist" contribution. 

To be specific, let us first consider the parton cross sec- 
tion for the 2 <Si/2 state. After some simple manipulation 
similar to Eq. jlj, we can reduce the amplitude into 
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bg{3 2 S h ) 



oc ^u h (p)-pT a v(p)e a (Z), (5) 
P-P 2 



where T a is the Zbb coupling, and e a (Z) is the polariza- 
tion vector of Z . Squaring the amplitude and summing 
over the color are standard. Less familiar is the sum over 
two polarizations of Uy 2 , yet can be accomplished with 
the aid of the following formula 



^2u f l(p)u'i(p)p fl p L , = ^(# + m)[(p-u) 



2 f\ 



(6) 



Note this resembles summing over polarizations of a 
Dirac spinor, up to a normalization. Finally we arrive 
at a succinct expression: 



t[bg(3 2 S 1/2 ) 



l§ira s {Mz) 1-47 



(1 - 2 7 )' 



T [bb], (7) 



where 7 = mj/M|, and Tq is the lowest order Z° partial 
decay width into bb: 



T [bb] = 



^^[^(1 + 27)+ fl »(l-47)], (8) 



where (3 = ^/l — 47, the vector coupling g v = — | + 
I sin 2 9\y and the axial- vector coupling g b A = — | . 

The HQR parton cross section for the A S 3 / 2 state can 
be derived in a similar way. It turns out t[bg(3 1 i S 3/2 )} = 

2T[bg(3, 2 Si/ 2 )), as expected from counting the polariza- 
tions. It is remarkable that they are simply To multiply- 
ing a constant. Of course, this simplification is mostly 
due to the vanishing of Fig.^i). 



The kinematic factor appearing in Eq. J7J suppresses 
the HQR contribution when the center-of-mass energy is 
near the bb threshold. Nevertheless, on the top of Z pole, 
setting this kinematic factor to 1 is an excellent approxi- 
mation for b and c. Since the HQR parton cross sections 
are finite when 7 = 0, it might be tempting to apply 
Eq. @ also to the light hadron production. However, 
the validity of HQR factorization in Eq. (JJJ is crucially 
based on the heavy quark dynamics, thus one should not 
erroneously use Eq. (JJJ) to describe Z° decay into light 
hadrons, by literally taking 7 — > 0. 

It may look surprising that the HQR parton cross sec- 
tions are 0(1) relative to r , not suppressed by m 2 /Mf 
as indicated in Eq. (JSJ). In fact, the explicit suppression 
from nib in Eq. (0 is compensated by p-v in Eq. ©. This 
non-suppression has important implication to the heavy 
flavor production in the hadron collision. By examining 
the bg recombination process gg — > bg(3) + b, one finds 
that the differential cross section da[bg(3))/dp 2 ~ 1/pj 
at p t 3> mi, a priori of the same order as fragmentation. 
In fact, it is equal to the differential cross section for gg 
annihilating into massless qq pair, up to a numerical con- 
stant. This finding diametrically challenges the general 
belief, that any non-fragmentation process at large p t will 
be suppressed by powers of 1/p 2 . 

To expedite the analysis, we choose the following 
linear combination of recombination factors: ^3 C ff = 
«3, 2 Si/ 2 )] + 2£[M3, 4 S 3 /2)]- We further define 6tot = 
J2~B^3cs[bg — ► B], where the sum is over all the lowest- 
lying b mesons and baryons. So the net contribution of 
this mechanism to the partial width of Z° to bb is 



AT[bb] = 32 f °'^ T [bb], 



(9) 



where we have doubled Eq. iQ to include the contribu- 
tion from CP conjugate of Eq. (JIJ. This is the crux of 
this work. From now on, we will use £3 as a shorthand 
for ^3tot for simplicity. 

The correction to the partial width of Z° into charm 
can be obtained by simply replacing £3 and To[&6] with 
£3 and Tofcc] in Eq. The heavy quark symmetry 

suggests that £3 = turrit, /m c s» 3 £3. However, a caveat 
is that, the symmetry-breaking effect may be large in 
the real world and practically we better treat them as 
two independent unknown parameters. 

The bg recombination can also make a correction to the 
hadronic width of W boson, through the decay channels 
W + — > bg + c(u) and W + — * eg + s(d) (Same reasoning 
can be also applied to t — > bg + W + ). Note it is not 
necessary for the recoiling quark to be heavy. The cal- 
culations are essentially the same as described above for 
Z° decay. However, unlike V z which is measured to per 
mille accuracy, is measured with a 3% error. There- 
fore the Z-pole observables are more sensitive probes to 
ascertain the effects of this mechanism. 
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TABLE I: Some selected electroweak variables, where R e = 
r ha d/r[e+ e -] and a° ad = 127rr[e+ e~]Y^ / {M%T%) . Listed 
are the latest measurements Q] and the SM predictions 0. 
Pull is defined as (O mcas - Ofit)/cr mcas . 





Meas. 


Predictions 


Pull 


r had [GeV] 


1.7444 ± 0.0020 


1.7429 ±0.0015 




T z [GeV] 


2.4952 ± 0.0023 


2.4972 ±0.0011 


-0.9 


^Ld [nb] 


41.541 ±0.037 


41.470 ±0.010 


1.9 


Re 


20.804 ± 0.050 


20.753 ± 0.012 


1.0 


Rt 


0.21644 ± 0.00065 


0.21572 ±0.00015 


1.1 


Rc 


0.1718 ±0.0031 


0.17231 ±0.00006 


-0.2 


T w [GeV] 


2.139 ±0.069 


2.0921 ± 0.0025 


0.7 



The discrepancy in A b FB is believed to either be an 
experimental problem, or originate from some tree-level 
new physics effects on Zbb vertex |2J ■ Our bg recombina- 
tion doesn't have much impact on the forward-backward 
asymmetries of b and c, since the role of this mechanism 
is simply rescaling both forward and backward cross sec- 
tions by a common factor, which cancels out in the ratio. 

From Table [Q we can see the general agreement be- 
tween the measurements and SM fits is rather good, 
therefore demanding £3 must be small. Naively, if we 
assume including the corrections in Eq. doesn't affect 
the values of SM fit parameters, then the relative vari- 
ations of the SM predictions listed in Table [I] are (For 
simplicity, we take £3 = 3 £3 temporarily) : 



gether with the SM parameters Mz, Mh, m t , a s (Mz) 
and Aa^ d . It will be very interesting to see how these 
parameters vary, and how much the quality of global elec- 
troweak fit may improve. 

To summarize, we have studied a new QCD mecha- 
nism that generates power corrections to partial widths 
of Z Q and W at 0(a s ), which have previously eluded 
the OPE analysis. The present precision of electroweak 
measurements requires these terms must be included. 

The phenomenological consequences of this mechanism 
to other collider experiments should be investigated. For 
example, we have noticed one unusual feature of bg re- 
combination in hadron collision - it is formally of the 
same order as fragmentation contribution at large pt- 
This is in sharp contrast to bq and bq recombination, 
whose effects are suppressed by Aqcd?™/?'* 0- I n addi- 
tion, because the dominant contribution to heavy flavor 
hadroproduction arises from the gluon-initiated subpro- 
cess, we expect bg recombination is more important than 
bq and bq recombination. 

Unfortunately, a rather (unnaturally) small £3 inferred 
from this work makes the color-triplet contribution virtu- 
ally invisible in hadron and ep collider experiments, due 
to much cruder consistency between data and the NLO 
QCD predictions. Since ^6 and £15 remain unconstrained, 
we hope they perhaps are much larger than £3 , so that the 
color-6 and 15 channels may result in noticeable effects. 

I thank C.-P. Yuan for discussing electroweak observ- 
ables. This work is supported by the National Science 
Foundation under Grant No. PHY-0100677. 



<5r had /r had w r(R b + 3R C ) 

STz/Tz = (Thud/Tz) ^had/That! 
°had — (1 — 2Fhad/rz) (JFhad/Thad 
SR e /R e = (ffhad/rhad 

SRb/Rt « r (1 - Rb - 3R C ) 
SR C /R C w t (3 - R b - 3R C ) 
ST W /T W w 3r/2(0.707GeV/r w ) 



S3' 

0.7 6 



3i 



-0.4 C; b 

?3' 

0.4 & 

hi 

0.7 & 



3- 



where r = 32ira s £3/9, and a s (Mz) = 0.12 is used. In the 
last row, we use T[W+ -> u.dj] « 0.707 GeV 0, 
and only include the Cabbibo-favored channels W + — > 
eg + s(d). To be compatible with most variables, £3 ~ 
10~ 3 looks reasonable. It will drive R e , Rb and Tw to- 
wards the correct direction, but deteriorate Tz, "had an< ^ 
R c . Note the correction to R c is much more significant 
than to Rb- A £3 as large as 0.01 becomes apparently 
unacceptable, since it would bring up the deviations in 
cr{J ad and R c to 6.4 a 01 and 1.8 a, respectively. 

Since many quantities are correlated in a complicated 
way, an unbiased strategy is to incorporate this mecha- 
nism into the global electroweak analysis. Consequently, 
the recombination factors £3 and £3, can be fitted to- 
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